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Abstract 

The crystal and magnetic structures of the specified CMR manganite system have been studied as 
a function of y = (0.2 — 1) across the metal-insulator (MI) transition, and of the oxygen mass ( 16 0, 
18 0). We quantitatively show how the polaronic narrowing of the carrier bandwidth and the crystal 
lattice micro-strains control the volume fractions of the mesoscopic ferro- and antiferromagnetic 
clusters. A well-defined dip in the transition temperatures and the suppression of all the types 
of long range ordering seen near the Mi-transition at y ~ 0.9 indicate a key role of the quenched 
disorder for the formation of the long-scale phase separated state. 

PACS numbers: 75.47.Gk, 61.12.Ld, 75.30.-m, 



1 



I. INTRODUCTION 



The presence of large 1 or giant oxygen isotope effects 2,3 in the colossal magnetoresistance 
(CMR) manganese oxides A!_ x A^Mn0 3 (A is rare earth element, A' is Ca, Sr, Ba) is now 
well established and it is proven to be an intrinsic property of these systems. Due to the 
strong electron-phonon coupling mediated by the Jahn- Teller effect the substitution of 16 O 
by 18 O can significantly change the effective exchange interaction between neighboring Mn 
ions 4 . In case of double exchange (DE) type of interaction the increase in the oxygen mass 
strongly reduces the effective charge carrier bandwidth, thus shifting the system towards 
insulating state, while the standard superexchange with large Hubbard energy gap is not 
supposed to be affected by the oxygen mass. Thus the value of the isotope effect on the 
metallic and insulating states of the above systems can reveal the underlying interatomic and 
electron lattice interactions responsible for particular magnetic, charge and orbital ordering, 
which are strongly coupled together because of the 3d-electron anisotropy and the mixed 
Mn 3+ /Mn 4+ valence. The low bandwidth manganite family (Lai_ y Pr y ) .7Cao.3Mn0 3 (LPCM 
hereafter) has the fixed optimal hole doping x = 0.3 and variable A-cation radius (r A ) 
that is linearly connected with the Pr-concentration y. The principal effect of decreasing 
(ta) is a decrease in Mn-O-Mn bond angle, which leads to the decrease in the electron 
transfer integral between the Mn-ions. The metal-insulator (M-I) boundary lays at the 
Pr-concentration between y = 0.86 and y = 1.0 5 ' 6 . The metallic state is conditioned by 
the well known double exchange and it is ferromagnetic (FMM), e.g. for y = 7 ' 8 . The 
magnetic structure of the insulating state can be both antiferromagnetic (AFMI) with the 
pseudo-CE type magnetic structure 9 ' 10 and ferromagnetic for the Pr-concentrations close to 
y — 1. The AFMI state is a Mott insulator state with the charge ordering of the Mn 3+ /Mn 4+ 
ions similar to the case of Lao.sCao.sMnOa, whereas the origin of the genuine ferromagnetic 
insulating phase (FMI) is not completely understood and it is characterized by a non- 
ordinary kind of orbital/charge ordering 11 ' 12 ' 13 ' 14 . The LPCM system can be converted from 
the metallic to insulating state by the isotope substitution 3 ' 15 . Similar giant isotope effect 
was observed in (La .5Ndo.5)o.67Cao.33Mn0 3 2 . The complete suppression of the ferromagnetic 
transition temperature cannot be attributed to a solely polaronic narrowing of the bandwidth 
within the DE model 1 ' 4 , but requires an inhomogeneous phase separated state, which was 
indeed directly observed by different experimental techniques in (La 1 _ y Pr y )o.7Cao.3Mn0 3 10 ' 16 
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and in the similar series La 5 / 8 _ y Pr y Ca3/ 8 Mn0 3 17,18 . The neutron diffraction study 15 has 
shown that the increase in the oxygen mass completely converts the dominating FMM 
metallic phase into the insulating charge ordered AFMI one in LPCM with y = 0.75. The 
presence of the long scale percolate phase separation in the manganites itself has been 
attracting special experimental and theoretical attention. There are two main concepts of the 
mesoscopically inhomogeneous state in the manganites. One standpoint is that the intrinsic 
quenched disorder enhances the fluctuations of the competing orders near the original bi- 
critical point 19 ' 20 ' 21 ' 22,23,24 . In another approach, the lattice distortions and the long-range 
strain 18 ' 25 ' 26 ' 27 are the dominant factors controlling the phase separation, similar to the one 
observed at the martensite type structural transition. Motivated by the spectacular giant 
isotope effect and the problem of phase separation we have undertaken a systematic neutron 
and synchrotron x-ray diffraction study of a series of the samples (Lai_ y Pr y ) .7Cao.3Mn03 
for y = 0.2 — 1.0 with 16 and 18 0. In the paper we present the detailed data on the 
structure and magnetic properties of the FMM, FMI and the charge ordered AFMI phases 
and discuss the key parameters conditioning the giant isotope effect and the phase separation 
in the LPCM-system. Earlier we have performed a diffraction study with another series of 
samples of the same system prepared by a different chemical route 10 , in which the effect 
of the oxygen isotope substitution has been studied only for one concentration y = 0.75. 
Moreover the different synthesis route used in the present study gives different microstructure 
parameters, which as we show below, have influence on the phase separation and hence on 
the oxygen substitution effect. 

The paper is organized as follows. In Sec. II the details of the sample synthesis, oxygen 
treatment and the experimental measurements are given. In Sec. Ill A the temperature 
dependences of the magnetic susceptibility xCO are analyzed. The ferromagnetic transition 
temperatures T c determined from x{T) are used as independent estimations of T c comple- 
mentary to ones determined from the neutron diffraction measurements. In Sec. Ill B the 
room temperature crystal structure and microstructure parameters as a function of y are 
compared for 16 and 18 0-samples. In Sec. Ill C the structure parameters are analyzed as a 
function of temperature between 15 K and 1200 K across the orbital order- disorder transition 
at ~ 700 K. In Sec. HID we show the evolution of the orbital and charge order as a function 
of doping y. In Sec. HIE the data on the temperature dependence of the magnetic state 
and the spatial phase separation are presented. The magnetic transition temperatures and 



3 



the magnetic ground state (the phase fractions, ordered magnetic moments) are obtained 
as a function of doping y and the oxygen mass. Sec. Ill F explains the connection between 
the phase separation and the lattice micro-strains (static variance of metric parameters). 
In Sec. Ill G the suppression of all types of ordering near critical concentration y c = 0.9 is 
discussed. This effect can be well interpreted in the theoretical model 19 ' 20 , which assumes 
the presence of compering states and quenched disorder near y c . 

II. SAMPLES. EXPERIMENTAL 

The crystal structure parameters of a LPCM sample for given y depend not only on y 
but also on the sample preparation procedure. For example, the spontaneous orthorhombic 
strain, thermal displacement parameters (reflecting also the local static disorder), the lattice 
micro-strains decrease as a function of annealing time and of the cooling time through the 
high temperature pseudocubic-orthorhombic structural transition. Thus, a precise study of 
a real LPCM system as a function of y and the oxygen mass can be done only within a series 
of samples prepared and thermally treated at the very identical conditions. The samples, 
which are denoted as "O-series" , were calcined at temperatures 1270-1570 K for lOOh from 
La20"3, Pr 6 On, MnO"2 oxides and CaC03 with 3 intermediate grindings. These samples we 
call "as-prepared". The final 18 0- and 16 0-samples were obtained via respective oxygen 
isotope exchange in closed quartz tubes in parallel under the controlled gas pressure slightly 
above 1 bar at 1270 K during 40h and then cooled down to 290 K with the rate 30 K/h. 
The 18 -samples had 80% of 18 0-isotope, measured by the weight gain after the oxygen 
exchange. The control weighing of the 16 0-sample gave the same mass within accuracy 
0.03%. The mass of each sample was about 2 g. The oxygen content in all the samples was 
determined by the thermogravimetric hydrogen reduction 28 and amounted to 3.003(5). The 
samples with y = 0.85 and 0.95 were prepared only with 16 0. One of the samples which 
has been studied earlier 10 ' 15 ' 29 , was also oxygen treated in the same way as all the samples 
from the O-series. This sample is from a series prepared by a synthesis from the aqueous 
solutions of the respective metals nitrates (denoted as "N-series"). 

The ac magnetic susceptibility xCO = x'(T) +ix"{T) was measured in zero external field 
with amplitude of the ac field 10 Oe and frequency 1 kHz using Quantum Design PPMS 
station. Neutron powder diffraction experiments were carried out at the SINQ spallation 
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source of Paul Scherrer Institute (Switzerland) using the high-resolution diffractometer for 
thermal neutrons HRPT 30 (A = 1.866, 1.494 A, high intensity mode Ad/d > 1.8 • 10~ 3 ), and 
the DMC diffractometer 31 situated at a supermirror coated guide for cold neutrons at SINQ 
(A = 2.56 A). All the temperature scans were carried out on heating. X-ray synchrotron 
diffraction at room temperature was done at the Material Sciences beam line (MS, SLS/PSI). 
The refinements of the crystal and magnetic structure parameters were carried out with 
FULLPROF 32 program, with the use of its internal tables for scattering lengths and magnetic 
form factors. 

III. RESULTS AND DISCUSSION 
A. Magnetic susceptibility 

The main use of the magnetic susceptibility data x{T) for the purposes of the present 
work is the determination of the ferromagnetic transition temperatures, complementary 
to the NPD data. The x(T) as a function of temperature is shown in Fig. 1 for several 
selected samples [since the x{T) are quite similar for the samples with close values of Pr- 
concentrations, not all the x(^)-curves are shown to avoid the plot overload]. This plot nicely 
illustrates the systematic decrease in the magnetic transition temperatures as a function of 
the oxygen mass. All the samples show (Fig. 2) the ferromagnetic type of the Curie- Weiss 
susceptibility x'(T) at high temperatures. The slope of the x f l (T) becomes non-linear 
below 250 K due to the presence of the antiferromagnetic correlations which are developed 
as the temperature decreases down to the Neel temperature T N ~ 150 K. The Curie- Weiss 
temperatures T C w were determined from the fit of the high temperature susceptibility to 
x'{T) = C/(T — T C w) m the temperature range (270 — 335 K) depending on y content. 
The Curie- Weiss temperatures amounted to Tew — 165 — 240 K for the 16 0-samples for 
y — 1 — 0.2. They are in agreement with the ferromagnetic transition temperatures Tq 
determined from NPD only for y = 0.2 and 0.5, while for the higher Pr-content Tew is 
significantly higher than T c . In the 18 0-substituted samples the T cw is decreased similar 
to the change in Tc, implying smaller effective electron transfer integral which mediates the 
ferromagnetic coupling for the heavier oxygen ion. The Curie- Weiss constant C = N A 2S(S+ 
1)a*b/3/c allows estimating the paramagnetic Mn-spin value S. The experimental values of 
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C amounted to 4.6 — 6 emu/mol depending on the y-content that correspond to the effective 
spin 2S ~ 5.2 — 6, instead of the expected value of 3.7 for the given Mn 3+ /Mn 4+ ratio. 
The enhanced value of the Curie- Weiss constant is a typical feature of the ferromagnetic 
manganite system and is attributed to the formation of the ferromagnetic clusters well above 
T<7 33 ' 34 . The quantitative analysis of C as a function of y is difficult because the Curie- Weiss 
constant value depends on the temperature interval used for the fit and the large variation 
of the To- The imaginary part of the susceptibility shows a maximum at a temperature 
very close to T c for all y (Fig. 1). The maximum is caused by the absorption due to the 
ferromagnetic domains formation with sufficiently short relaxation times r ~ l/(27r/) close 
to Tq, where / is the external field frequency. This maximum is asymmetric and can be 
very broad, so the temperature of the maximum cannot be used for quantitative estimation 
of the Tq. A good way of the evaluation of the ferromagnetic transition temperature is the 
analysis of the derivative of the magnetic susceptibility dx'/dT(T) which has a well defined 
minimum as shown in the insert of Fig. 2. The temperature of the minimum T m j„ will be 
taken as an independent, complementary to the NPD data, estimation of Tc- 

B. Crystal structure and microstructure 

The crystal structure for all the compositions y = 0.2 — 1.0 at all temperatures is well 
refined in single phase in the orthorhombic space group Prima with the standard for these 
compounds structure model (see e.g. 29 ). In general, the low temperature state is spatially 
mesoscopically phase separated, but the lattice constants of the phases are too close to each 
other and the two phase refinement is not reliable. However, there are cases where the 
structure separation is explicitly visible, e.g. in the system Ro.sSro.sMnOs (R = Sm, Td, Tb) 
the metrics of the phases are sufficiently different and both phases can be identified in the 
neutron diffraction experiment 35 . The examples of the x-ray synchrotron diffraction pattern 
and the neutron diffraction pattern and their refinements are shown in Fig. 3 and Fig. 4, 
respectively. All the room temperature (T = 290 K) crystal structure parameters except for 
the thermal displacement parameter do not depend on the oxygen mass within the accuracy 
0.2%. The lattice constants, which are determined with the best accuracy (Fig. 5(a)) are 
slightly systematically smaller for the 18 0-samples: for y < 0.5 the unit cell volume V is 0.2% 
smaller, for the y > 0.5 the volumes V match within 0.02%. The only structure parameter 
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which is noticeably and systematically changed by the oxygen isotope substitution is the 
thermal displacement parameter B. The oxygen thermal parameter is slightly decreased, 
while the Mn and A-cation ones are pronouncedly increased as a function of oxygen mass 
(Fig. 5(b,c)). These changes agree with the expected behavior of the B ~ (A;^T)a; _2 /M, 
where M is the atom mass, the u~ 2 is averaged over the phonon density of states Z(oo) 3e . 
The parameter u~ 2 should increase as a function of an average atom mass provided that the 
elastic force constants are not changed. Hence the thermal parameters of Mn and A-cation 
should increase, while the oxygen 5-value behavior cannot be simply predicted because 
both uo~ 2 and M are changing. Figure 6 illustrates systematic distortion of the structure 
as a function of the Pr-concentration y. The Mn-O-Mn bond angles are systematically 
decreased conditioning the change from the metallic to insulating state at the concentration 
above ~ 0.85. At the same time the orthorhombic strain r is substantially increased with 
the increase in y. The values of strain r for the "as-prepared" samples (also shown in the 
Fig. 6a) are significantly larger stressing the point that the systematic study of the present 
system as a function of y and the oxygen mass is only possible within a series of samples 
prepared and treated in identical way. 

The microstructure of the samples was obtained from the analysis of the Bragg peak 
broadening (anisotropic in general) as a function of scattering angle (26*) by the following 
procedure. The Bragg peaks were described by the pseudo-Voight distribution function 
with the Gaussian (r G ) and Lorentzian (r L ) components. The total peak width T was 
calculated from the individual components using a 5th order polynomial approximation as 
described in Ref. 37 and implemented in the FULLPROF code 32 . The angular dependence of 
the peak width was approximated by the Cagliotti-type function Tg(26) = ((£/ tan 2 9 + 
a%) + VtanO + W + o-ig/cos 2 ^ 1 / 2 ) and T L (26) = (Xtantf + a L ) + (Y + a 1L )/cos6 for 
the Gaussian and Lorentzian components respectively. All the widths correspond to the 
full width at half-maximum (FWHM) and the angles are measured in radians. U, W , U 
and Y were fixed to the values given by the instrumental resolution function. aa(26) and 
<7£,(20) are the Rietveld refined Gaussian and Lorentzian components related to the sample 
micro-strains. G\q and o\l are the refined components related to the finite size broadening. 
The total widths both for the micro-strain (o (29)) and the finite size (a±) broadening were 
calculated from the individual Gaussian and Lorentzian refined components by means the 
above mentioned 5th-order polynomial approximation. The apparent sizes (L) are connected 
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with the broadening a\ by well-known Debay formula a± = X/L. The micro-strains or 
the static fluctuations of crystal lattice constants are determined from the variance of the 
metric parameters (A,B,C) of the reciprocal space. The square of the reciprocal distance 
1/d for a Bragg reflection with Miller indices (hkl) for the orthorhombic singony reads 
1/d 2 = M hkl = Ah 2 + Bk 2 + CI 2 . The variance of M hkl is related to the widths of the 
diffraction peaks as a (29) = a(Mhki)/Mhkit&n9. General description of the anisotropic 
peak broadening is given in Ref. 38 . In our case we consider the isotropic micro-strains 
(5d st /d), where d stands for a, b or c lattice constants, and an additional micro-strain 
along a- axis (Sd st /d)i W0 ] that have the contributions to the reciprocal distance variance 
<j(Mhki) = 2(5 d st I ' d) M^i and 2(5d st /d)[ 100 ]Ah 2 , respectively. The presence of the anisotropic 
broadening is nicely illustrated by the two doublets shown in the insert of Fig. 3. 

The sample related peak broadening is comparable with the resolution of neutron diffrac- 
tion instrument HRPT, therefore we have performed a comparative ultra-high resolution 
synchrotron x-ray diffraction study of several samples at room temperature. The examples 
of the reduced sample related peak widths 5d/d = r samp i c / tan(0) as a function of scattering 
vector q = 2ir/d are shown in the insert of Fig. 4. The width r samp i e includes the isotropic 
micro-strains and the finite size broadening components. One can see that the experimental 
broadening is dominated by the micro-strain effect, which gives a constant term to 5d/d(q), 
whereas the finite size effect gives a oc l/q term. The apparent sizes L are larger than 2-10 3 A 
giving a relatively small contribution to the broadening (assuming 1/L — increases the 
refined micro-strain values by less than 10%). The micro-strains refined from the neutron 
and synchrotron data agree within the statistical error bars. This implies that in spite of 
lower resolution of the neutron diffraction data (Sd/d ~ 2 • 10~ 4 and > 1.8 ■ 10~ 3 for the 
synchrotron x-ray and neutrons respectively) we can reliably extract the microstructure pa- 
rameters from the NPD data. The lattice micro-strains allow us to characterize the amount 
of quenched disorder and the presence of the charge ordering as discussed in the following 
sections. 

There are two transitions in (Lai_ y Pr y ) .7Cao.3Mn03 associated with the changes in the 
crystal structure. The first transition at high temperature is a fist order isostructural phase 
transition from the pseudo-cubic to the orthorhombic orbital ordered structure similar to the 
one observed in the stoichiometric LaMn0 3 39 . The second transition at lower temperature is 
related to the charge ordering of the Mn-ions and is accompanied by the lowering symmetry 
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to the monoclinic one and the unit cell doubling along the orthorhombic a-axis. We have 
studied the crystal structure of several compositions of LPCM up to the high temperatures 
across the transition to the pseudocubic phase. 

C. Orbital order-disorder transition 

The study the crystal structure of LPCM up to the transition to the pseudocubic phase 
allowed us to get a wide range of the changes in the crystal structure, which helped to 
understand the temperature behavior of the crystal structure at low temperature. Figure 7 
shows the crystal lattice constants and the isotropic micro-strain parameter (5d st /d) in the 
LPCM sample with y=0.75. The latent heat at the transition temperature Tj-t (insert in 
Fig. 7) points to the first order phase transition. Above the transition the crystal structure 
is metrically cubic and the refinement in the profile matching mode in the cubic metrics 
gives the same quality of the fit (x 2 ) with the lattice constant shown by rhombs in Fig. 7a. 
In addition, the Mn0 6 octahedron gets regular above the transition. The structure is thus 
becoming pseudo-cubic. This has inspired us to search for a possible alternative cubic crystal 
structure description, similar to the attempt undertaken in 39 . The diffraction pattern at 800 
K is indexed in the cubic cell with a = 7.76 A. The list of possible space groups with their 
numbers in brackets, which does not contradict the set of observed diffraction peaks contains 
seven possibilities: P2 : 3 (198), P4 2 3m (208), P23 (195), Pm3(200), P432 (207), P432 (215) 
and PmZm (221). Our trial structure solutions (carried out with the program FOX 40 ) did 
however fail to find a suitable candidate structure model. As it was actually pointed out 
in 39 , the failure could have been predicted by the even earlier studies by Glazer 41 , who did 
not deduce any primitive cubic space group for the structure with lattice constant 2a p for 
the distorted perovskite which could have originated from the coherent octahedra tilts. 

The micro-strain value (5d st /d) is markedly increased below the transition temperature 
(Fig. 7b). The lattice constant mismatch in the adjacent crystal domains and/or twins 
naturally creates the internal micro-strains below the transition Tj_ T . This is important to 
note because the phase separation at low temperature depends on the value of the micro- 
strains, which will be inevitably present in the LPCM system due to this structure transition. 

Figure 8 shows the bond lengths between manganese and oxygen atoms in the MnOg- 
octahedra as a function of temperature. Mn-Ol is the bond directed to the oxygen in the 
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position shifted along y-direction with respect to the Mn-atom. Mn-021 and Mn-022 are 
the bonds to the oxygen atoms located roughly in the (ac) plane. The transition at Xj_ T is a 
Jahn- Teller type of transition connected with the lifting up of the e g orbital degeneracy. The 
filling of the e g orbitals \e g ) = cos(#/2) \z 2 ) + sin(#/2) \ x<2 ~ U 2 ) can be determined from the 
value of an angle 9, which is the polar angle in the plane of normal coordinates (Q2, Qz) 42 - 
The coordinates are given by the distortion of the MnC>6 octahedron Q3 = (l/v / 2)(/ — s), 
Q2 = (2/v / 6)(2m — I — s), where s, m, I are the short, medium and long Mn-0 bond lengths, 
respectively. One can see from the Fig. 8 that above Tj_ T all three bonds are equal, implying 
that both e 3 -orbitals are equally populated. Below the transition two of them remain equal 
(m = s), that gives the angle 9 = 27r/3. This angle corresponds to the complete filling of the 
y 2 -orbital or for another choice of axes to equivalent z 2 - or x 2 - orbitals [they are related to 
each other by 2n/3 rotation in the (Q 3 , ) -plane] . In the Pnma structure it corresponds to 
the antiferrodistorsive type of orbital ordering (00) schematically shown in the Fig. 8. It is 
interesting to compare with the case of LaMnOs, where below the transition the z 2 -orbital 
is ~ 65% filled, while above the transition the angle 9 amounted to 9 ~ n/2 (Q 3 = 0), 
which also corresponds to the equal filling of the e g orbitals 39 . In the later case the Q2-type 
of the octahedron distortions remain in the high-temperature phase due to the interactions 
between the distortions. In our case the interaction between the Jahn- Teller active Mn 3+ 
ions is weakened due to the dilution by 30% of the statistically distributed Mn 4+ ions and 
thus the MnOe octahedron behaves like a free molecule with equally populated e 9 -orbitals. 
The difference in the in-plane bond lengths Mn-02 in the 00 state at room temperature is 
an order of magnitude smaller in the case of LPCM (cf. 0.27 A vs. 0.02 A). This abrupt 
decrease in the bond length mismatch as a function of the Mn 3+ fraction is an additional 
evidence of the cooperative nature of the Jahn- Teller effect. In the Lao.sCao.sMnOa where 
the ratio of Mn 3+ : Mn 4+ is "50:50" the Mn0 6 -octahedron remains almost regular at room 
temperature 43 . 

D. Charge order transition 

The LPCM system undergoes a transition to the charge ordered state (CO) at Tco > 
TV 9 ' 15 . The transition is revealed by the appearance of weak superstructure peaks reflecting 
the doubling of the unit cell along a-axis and the lowering of the symmetry probably to 
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the monoclinic P2 1 /m one similar to the case of Lao.5Ca .5Mn0 3 43 . However, the deviations 
from the average Prima structure are small and do not allow a reliable detection of the charge 
ordering. Better structure indicators of the CO-transition are the increase in the Mn-0 bond 
lengths mismatch and the appearance of the anisotropic lattice micro-strains in the CO state. 
One can see (Fig. 8) that below ~ 250 K the two equal Mn-0 bond length become different. 
We believe that this is a manifestation of the charge ordering which is schematically shown 
as insert in Fig. 8. This type of Mn 3+ /Mn 4+ and orbital ordering in the (ac) plane is dictated 
by the Goodenogh-Kanamori rules and the low temperature magnetic structure (PCE). In 
Ref. 9 there was proposed that in Pr .7Cao.3Mn0 3 the sites shown with the z 2 -type orbital 
are fully occupied by Mn 3+ , while those shown by circles are occupied by Mn 3+ : Mn 4+ with 
the ratio 40 : 60. Since the structure is described in the average Prima space group the 
atom positions are averaged over the symmetry elements of Prima. The orbital ordering in 
the charge ordered phase would give roughly equal average in-plane Mn-02 bond length. 
However, even in an ideal case of CO in Lao.5Ca .5Mn0 3 the in-plane Mn-02 bond lengths 
become different in the CO-phase at low temperature in the Pnma-description 43 . The in- 
plane Mn-02 bond length mismatch in the LPCM is larger apparently due to the larger 
Mn 3+ fraction. An evident indicator of the charge ordering is the mismatch of the short 
Mn-022 and Mn-Ol bond lengths. In the antiferrodistorsive 00 state all the z 2 -orbitals 
are in the ac-plane, while below the CO-transition the z 2 -orbitals located at "Mn 4+ " sites 
should be oriented perpendicular to the (ac)-plane to stabilize the ferromagnetic coupling 
between the neighboring planes along 5-axis of the PCE-type antiferromagnetic structure. 
This is reflected as a decrease in the short Mn-022 with respect to the Mn-Ol bond below 
~ 250 K (Fig. 8). Figure 9 shows the temperature dependences of the bond lengths as a 
function of Pr-concentration y. One can see that the difference of the Mn-Ol and Mn-022 
bond lengths is drastically decreased near y = 0.9, reflecting the suppression of the charge 
ordering. In the metallic predominantly ferromagnetic state (y < 0.5) CO is suppressed as 
well. Another parameter which can be considered as a structure indicator of the CO state is 
the anisotropic micro-strain (Sd s t/d)[i 00 ] along a-axis, along which the doubling of the unit 
cell occurs in the CO phase. This implies that the CO/OO is not perfectly homogeneous 
as in purely "50:50" composition Lao.5Ca .5Mn0 3 but there is a distribution of domains 
with different extent of CO thus leading to the dispersion of the lattice constant a. The 
anisotropic micro-strain is pronouncedly increased on cooling for y > 0.7, but has a nonzero 
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value at room temperature 290 K as shown in Fig. 10. A remarkable feature of (Sd st /d)[ W0 ] 
is a local minimum around y c = 0.9 pointing to a suppression of the CO concomitantly with 
the suppression of the AFM and FM ordering as discussed below. The isotropic micro-strain 
value (5d st /d)i so (also shown in Fig. 10) doesn't exhibit any peculiarity at y c ruling out any 
kind of sample effect. 

E. Magnetic state 

The temperature dependences of the magnetic Bragg peak intensities I(T) allow accurate 
determining of the transition temperatures to the AFM and FM states (T N , T c ). A typical 
diffraction pattern is shown in Fig. 4, where the tick rows indicate the nuclear and mag- 
netic phases. The magnetic structure consists of AFM pseudo-CE component with Mn-spin 
directed along b and FM component in (ac)-plane similar to the one reported in 9,10 . In 
addition, the Pr-spins also get ferromagnetically ordered parallel to the FM-component of 
Mn-spins. The refined low temperature magnetic moments are presented in Table I. Figure 
11 shows several examples of the temperature dependences of the selected magnetic Bragg 
peaks I(T). The 16 and 18 samples for each ^-concentration were prepared in the same 
amounts and measured in the identical experimental conditions and thus the integrated in- 
tensities shown in Fig. 11 were normalized only to the neutron monitor (number of incoming 
neutrons). Above 100K the samples with y > 0.7 are in purely AFM state that is not af- 
fected by the oxygen substitution. The Neel temperatures T N are the same and the whole 
sample volume is occupied by the AFM phase, as illustrated by Fig. lib. The ferromagnetic 
phase is suppressed by the increase in the oxygen mass that is reflected by the splitting of 
the AFM integrated intensities for the 18 and 16 samples roughly below 100K. 

One sees that the AFM transition is very well defined by the abrupt vanishing of the 
respective Bragg peaks. The ferromagnetic state is destroyed in a more smoother way. 
Some fraction of the sample volume remains ferromagnetic above the transition as one can 
see from the smeared shape of the FM Bragg peak intensities I(T) near T c . To determine the 
transition temperatures in an unbiased systematic way we performed least-square fits of the 
experimental I(T) in assumption of the Gaussian distribution of the transition temperatures 
and a phenomenological power-low function for the spontaneous magnetization: 
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I(T) = h r - 7 2— e (-^) 2 « (1 _ (T/rmi - H(T - r)) rfr + J , (1) 
Jo \i2ixbT c 

where if is a step-like Heaviside function, Tc and 5Tq are the mean transition temperature 
and r.m.s. variance, respectively. I\ are I are the refinable parameters. For the compositions 
with y > 0.5 the parameters a and (5 were fixed as indicated in Table I. Otherwise, the fit 
parameters got unreasonable values. To reduce a number of the refined parameters the 
experimental I(T) curves for 18 O and 16 0-samples were processed in one combined fit with 
the same parameters a, (3 and I . The results of the fits of the ferromagnetic /(/) to the 
formula (1) are shown by thick lines in Fig. 11. For comparison, the thin lines for the sample 
with y = 0.5 show the calculated I(T) curves deconvoluted from the Gaussian smearing. The 
magnetic transition temperatures are summarized in the Table I and in Fig. 12. The Curie 
temperatures Tc are in reasonable agreement with the ones determined from the magnetic 
susceptibility data that are also given in the Table I. One can see from the Fig. 12 that 
the isotope substitution has a large effect on the T c , whereas the Neel temperature T N is 
invariant as one would expect for the double exchange type of FM and the superexchange 
AFM interactions, respectively. 

The AFM phase has been detected at all concentrations, including y = 0.2. T N for 
the sample with y = 0.2 was not determined, because the AF-fraction is very small [the 
effective AFM moment at T = 15 K is m A = 0.24(6)// B ] and the temperature scan would 
require extremely long measurements. For y > 0.7, the AF ordering precedes appearance of 
the FM phase with the temperature lowering. The FM phase appears spatially separated, 
reducing the sample volume occupied by the AFM phase and thus reducing the slope of the 
Iaf(T) roughly below Tc- The sizes of the FM- and AFM-do mains are mesoscopically large 
(> 10 3 A) estimated in the similar way as described in Ref. 10 . The same effect of phase 
separation is present for y = 0.5, where T c > T N (Fig 11a) implying that there are two 
characteristic energy scales connected with the ferromagnetism. One (T c ) is given by the 
double exchange mechanism with the strong dependence on y. The second characteristic 
temperature is about 130 K where the phase redistribution in favor of the FM state is 
triggered on. 

Figure 12 shows the effective FM m F (y) and AFM rriA{y) magnetic moments of Mn- 
ions at base temperature (10-15 K) as a function of y. is the average value of the 
moments for two sublattices ttiai and presented in Table I. The effective moment is 
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proportional to the product of the real moment and square root of the volume occupied by 
the ordered phase. Since the crystal structure is well refined in a single phase approximation 
one cannot determine from a single diffraction measurement whether the magnetic state 
is a homogeneous one. However, it becomes possible from the analysis of the values of 
the effective magnetic moments for all the concentrations y for both 16 O and 18 0-samples. 
Figure 13 shows m\{m 2 F ) for all the samples with the Pr-concentration y as a parameter. 
The line is a least-square fit to the linear function m 2 A = M\(l — m 2 F /M F ) for all y except 
y = 0.9 and y — 1 with 18 0. This linear relation is virtually possible only in assumption 
of the phase separated state. Otherwise, one would have to suppose quite unreasonable 
sizes of the magnetic moment as a function of y and the oxygen mass. This linear function 
assumes that the magnetic moments in the AFM and FM phases keep the constant values 
Ma = 2.26(l)/iB, Mp = 3.57(2)//b over the whole phase diagram, but the phase volumes 
are synchronously redistributed as indicated by top-x and right-y axes of Fig. 13. The size 
of the FM moment well corresponds to the average spin value of Mn 3+ /Mn 4+ ions 3.7 
whereas the AFM moment is substantially lower. It can be connected with a "defect" charge 
ordering that causes a frustration in the magnetic exchange between the neighboring (ac) 
planes. With the help of the Fig. 13 one can disentangle the effect of the FM/AFM volume 
redistribution (the points slide along the straight line) and the effect of the suppression of the 
magnetic ordering (the points lay below the line). The increase in the oxygen mass simply 
shifts the balance between the phases towards the AFM one for the "metallic" part of the 
phase diagram (y < 0.8), similar as the increase in y does, implying that the effect of the 
isotope substitution is the polaronic narrowing of the bandwidth. The FM phase is metallic 
and if its volume decreases below the percolation limit ~ 16% 4 then the sample should 
become an insulator. In the "insulating" part of the diagram (y > 0.9) the system behavior 
is very different. Firstly, the ferromagnetic phase is not any more metallic because its 
fraction, e.g. in the sample with y — 1 is far above a percolation threshold corroborating the 
insulating character of the FMI phase. Secondly, the increase in the oxygen mass suppresses 
only the FMI phase, both Tc and the effective magnetic moment, leaving the AFM phase 
fraction intact. One can speculate that the FMI cannot be converted to the AFM phase 
because it is already insulating and thus the further decrease in the electron transfer integral 
doesn't necessarily shift the system towards the AFM state. We suppose that those spins that 
stop contributing to the FM phase become disordered thus reducing the effective moment. 
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However we cannot distinguish whether a spatially separated disordered phase is formed or 
some spins inside the ferromagnetic host get frustrated. The presence of the isotope effect 
on the FMI is remarkable itself, because it rules out an ordinary superexchange mechanism 
of FMI. This has been noted in Ref. 11 where the isotope effect on Tq was found in insulating 
Pr 1 _ x Ca x Mn0 3 , and the genuine FMI state was proposed to appear due to some CO/OO 
superstructure or density wave giving an energy gap which is much smaller than the Hubbard 
on-cite repulsion. 

F. Effect of micro-strains on phase separation 

The specific low temperature FM/AFM volume fractions shown in Fig. 13 do not repre- 
sent themselves universal for LPCM system values that depend only on y. Earlier we studied 
virtually the same series of samples only with 16 0(N-series) 10 , but prepared by the different 
chemical route than the one explained in section II. The phase redistribution process in 
the N-series samples (Fig. 8 of Ref. 10 ) triggered at about T c on cooling goes much "easier" 
leaving predominantly (> 90%) FM ground state for y > 0.8, while in the O-series samples 
the FM fraction is decreased to 73% for y = 0.8. The increase in the oxygen mass leads to a 
complete phase volume redistribution in favor of the AFM phase in the N-series sample with 
y = 0.75 15 . In the case of the O-series samples, the phase balance is also definitely shifted 
towards the AFM state (Fig. 13), but none of the 18 0-sample shows purely AFM ground 
state. We have repeated the oxygen treatment and oxygen isotope exchange with the N- 
series sample in the same way as we did with the samples of the present O-series, including 
the back exchange and got essentially the same results as in 3,15 . It is believed that the qual- 
itative difference between the two series originates from smaller grain sizes of the starting 
reagents in case of the N-series that could give different micro-structure of the sample crys- 
tallites. We have undertaken a quantitative micro-structure comparison of two exemplary 
samples with y = 0.75 from each series also using ultra high resolution synchrotron x-ray 
diffraction (as described in sec. Ill B). Both samples are well refined with the close values 
of structure parameters. The noticeable difference between the samples from the different 
series is their micro-structure. The insert of Fig. 4 shows the g-dependence of the reduced 
peak broadening 5d/d in the two samples from both series. One can see that the broadening 
is larger in the O-series sample. The main contribution to 5d/d comes from the micro-strain 
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effect amounted to Sd st /d = 0.16% and 0.23% for the N- and O-series sample respectively. 
Inside the O-series the isotropic micro-strain varies only slightly as shown in the Fig. 10. 
Thus, the key difference between the two series of the LPCM-samples conditioning their 
different ground states originates from the different values of (5d st /d) motivating to propose 
that the phase separation is favored by the micro-strains. The value of the micro-strains 
can be considered as a quantitative measure of the disorder in the crystal lattice. Many 
studies of the epitaxial films of manganites have demonstrated that the relative fractions of 
the FM/AFM phases can be strongly modified by both morphology and strains of the films 
(e.g. 44 ). The strong micro-strain dependence of the phase separated volume fractions in our 
case suggests that the value of the micro-strains (5d st /d) in the bulk material should also be 
considered as an additional parameter governing the formation of the phase separated state. 

G. Suppression of ordering near y c = 0.9 

The Curie temperature T c is decreased as a function of Pr-concentration y due to the 
smooth decrease in the Mn-O-Mn bond angle (Fig. 6b), but near the critical concentration 
y c = 0.9 one can see a local minimum (Fig. 12c), which is unexpected from the behavior of 
crystal structure parameters. The Neel temperature has a pronounced minimum at y c , as 
well. The total volume fraction occupied by both FM and AFM ordered phases decreases 
by 25% at y c in comparison with the other samples (Fig. 13). The charge ordering that is 
concomitant to the antiferromagnetic PCE-type ordering is also suppressed near the critical 
concentration y c (section HID). 

The suppression of the magnetic transition temperatures and the ordered phase volumes 
at y c = 0.9 strongly resembles the phase diagram calculated in Ref. 19 . The authors of 
Ref. 19,20 have shown that introducing of the quenched disorder into an Ising system with 
two competing ordered states gives an effect of a suppression of the ordering temperatures 
and mesoscopic phase separation in the vicinity of the bi-critical point. The experimental 
situation in our case is more complicated. Instead of two phases there are three distinct 
phases: AFMI phase and concomitant CO phase, FMM metallic phase for y < y c and FMI 
insulating phase for y > y c . The phases, which would correspond to the phases separated by 
the critical point are the FMM metallic from one side and both FMI and AFM insulating 
phases from another side of MI transition. In the spirit of the above model one would 
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expect that there is also a phase separation into the FMM and FMI phases near y c , which 
are indistinguishable from the present experimental data. The quenched disorder in the 
LPCM-system is naturally present due the dispersion of the A-cation radius amounting to 
a a ~ 0.01 A near y c (this value would nominally correspond to quite large Pr-concentration 
spread Ay ~ 0.4). However, since the radii of the Pr and Ca ions are almost the same 
the dispersion a a limits to zero for y — 1, while the phase separated state is still present. 
The role of the quenched disorder can play the Jahn- Teller strain from Mn 3+ centres as 
proposed in Ref. 45 . In addition, the quenched disorder can originate from the chemically 
inhomogeneous distribution of the A-cations given by the sample preparation procedure. 

We cannot exclude also that the lattice distortions and the long-range strain 18 ' 25 can be 
the reason for the phase separated state. We note that the suppression effect was not found 
in a detailed study of similar system Pr .65(Ca y Sri_ y )o.35Mn03 22 , where all transition tem- 
peratures were continuous through the critical point. The later system has almost the same 
average A-cation radius but much larger (~ 4 times) variance (Ja, which could be a possi- 
ble reason of the different behavior. We suppose that to compare the different manganite 
systems with respect to the effect of the suppression of the ordering near y c one should also 
compare the extent of quenched disorder, which can be quantitatively characterized by the 
value of the isotropic micro-strains. 

IV. CONCLUSIONS 

The magnetic and crystal structures have been studied in the series of samples 
(Lai_yPr y )o.7Ca .3Mn03 as a function of the Pr-concentration (y = 0.2, 0.5, 0.7, 0.75, 0.8, 
0.85, 0.9, 0.95, 1.0) across the metal-insulator (MI) transition around y c = 0.9 and as a func- 
tion of the oxygen mass ( 16 0, 18 0). The ground magnetic state is mesoscopically separated 
into antiferromagnetic insulating (AFM) and ferromagnetic (FM) phases for all values of y. 
The ordered magnetic moments of the Mn-ions keep the constant values M AFM = 2.26(l)/x B , 
M FM = 3.57(2)/i B over the whole phase diagram, whereas the phase fractions are changing 
as a function of y. The conductivity in the "metallic" part of the phase diagram (y < y c ) is 
metallic because the FM phase is metallic and its fraction is above the percolation threshold. 
In the "insulating" domain (y > y c ) the fraction of the FM phase is significantly above the 
percolation threshold (e.g. for y = 1.0 it amounts to 46%) implying that the FM phase is 
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insulating for y > y c . 

The increase in the oxygen mass has effect on both the values of T c and the phase 
fractions. The Curie temperatures of both metallic and insulating FM phases are decreasing 
in a similar way suggesting that the origin of the ferromagnetism in the insulating state 
should be connected with a slow electron hopping process similar to the double exchange 
as discussed in 11 . The Neel temperature is not affected by the oxygen mass as expected for 
the AFM superexchange. The increase in the oxygen mass always suppress the FM phase 
fraction. However, the effect is different for the metallic and insulating FM phase. For 
y < y c , the oxygen substitution acts similar to the increase in y, namely the FM phase 
volume fraction is decreased by the same amount as the AFM volume fraction is increased. 
For y > y c the ferromagnetic insulating phase is suppressed but the AFM phase is not 
changed at all. Thus, the narrowing of the bandwidth due to the decrease in the A-cation 
radius and the polaronic narrowing act very differently in the "insulating" part of the phase 
diagram. 

The ratio between the AFM and FM phase fractions is not uniquely a function of the 
effective bandwidth given by y and the oxygen mass, but is also affected by the value of 
internal micro-strains that favor the formation of the phase separated state. We propose 
that the value of the micro-strains can be used as a universal quantity characterizing the 
extent of the quenched disorder, which can have different origins in the manganites. 

The phase diagram of (Lai_ y Pr y )o.7Cao.3Mn0 3 has a critical concentration at y c = 0.9, 
where both Neel T N and Curie T c transition temperatures, the ordered phase volumes as well 
as the charge ordering concomitant to the AFM phase are suppressed. This finding together 
with the presence of the phase separated state strongly support the theoretical study of the 
competition between two ordered states in the presence of quenched disorder 19 ' 20 , implying 
that the quenched disorder plays a key role in the formation of the phase inhomogeneous 
state in manganites. 
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TABLE I: The effective ferro- and antiferromagnetic moments at T < 15 K and the Neel T/v and 
Curie Tq transition temperatures in the (Lai_ y Pr y )o.7Cao.3Mn03. Tq were determined from the 
fits of neutron diffraction intensities to formula (1) (ND) and by the inflection point of magnetic 
susceptibility (%). 5Tc is r.m.s. variance of the Gaussian distribution of the transition temper- 
atures. See the text for details. The values of the moments are given in ^b- The Mn-magnetic 
moments of the pseudo-CE AFM structure are thai and mA2 for the propagation vectors k\ = [0 
1/2] and A;2=[l/2 1/2], respectively, mp is the FM moment of Mn-ion. mp r shows the magnetic 
moment of Pr-ion. 
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T(K) 

FIG. 1: (color online) Real x' an d imaginary x" parts of the ac magnetic susceptibility are shown 
as a function of temperature for several samples with the Pr-content y indicated in the plots. 16 
and 18 0-samples are shown by solid and dashed lines respectively, x" f° r V = 0-5 (which is similar 
to one for y = 0.2) is not shown to avoid plot overload. 
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FIG. 2: (color online) x' -1 as a function of temperature for the Pr-content y indicated in the plots. 
The dotted lines show the fits of high temperature part of x' to the Curie- Weiss low. The refined 
Curie temperatures Tew are indicated. The insert shows two examples of the first derivative 
dx'/dT as a function of temperature. 16 and 18 0-samples are shown by solid and dashed lines, 
respectively except for the y=1.0 ( 16 0) sample, which is shown by circles. The solid line in the 
insert for this sample shows example of the parabolic fit. The refined temperature of the minimum 
is indicated in the insert. 
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FIG. 3: An example of the Rietveld refinement pattern and difference plot of x-ray synchrotron 
diffraction data (MSP/SLS, PSI, A = 0.92 A) for the sample with y=I.O at room temperature. 
The insert illustrates anisotropic line broadening along a-axis. 
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FIG. 4: An example of the Rietveld refinement pattern and difference plot of the neutron diffraction 
data for the sample with y = 0.9 at T=15K. The most intensive magnetic peaks are indicated by 
Miller indices. The insert illustrates the reduced FWHM of the Bragg peak broadening (Sd/d) as a 
function of q for two samples with y = 0.75 from O- and N-series at room temperature measured by 
neutron and x-ray synchrotron diffraction. (Sd/d) was calculated from the refined FWHM-width of 
the Bragg peaks after deconvoluted with the instrument resolution function. 5d s t/d are the refined 
isotropic micro-strain value. See the text for details. 
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FIG. 5: (color online) (a) The lattice constants (sp. gr. Pnma) and (b,c) isotropic thermal 
parameters at room temperature T = 290 K. 16 O and 18 0-samples are shown by closed and open 
symbols respectively. 
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FIG. 6: (color online) (a) The orthorhombic strain r and (b) Mn-O-Mn bond angles as a function 
of Pr concentration (sp. gr. Prima) at room temperature T = 290 K. 16 and 18 0-samples are 
shown by closed and open symbols respectively. The triangles in (a) show the orthorhombic strain 
in the "as prepared" samples for comparison. The lines are guides for the eyes. 
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FIG. 7: (color online) (a) Lattice constants and (b) isotropic micro-strains as a function of tem- 
perature in the (Lai_ y Pr y )o.7Cao.3Mn03 (y=0.75). The insert shows the differential scanning 
calorimetry signal with the same x-axis. The minimum indicate the presence of the latent heat of 
the first order phase transition at the temperature Tj_t- 
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FIG. 8: (color online) Bond length Mn-Ol and Mn-02 being approximately along and perpendic- 
ular 6-axis, respectively as a function of temperature in (Lai_ y Pr y )o.7Cao.3Mn03 for y=0.7. The 
structure undergoes a change from the orbitally disordered pseudocubic phase to the antiferrodis- 
torsive type of orbital ordering (00) below 800K and to the charge ordered (CO) state at lower 
temperatures. The sketches of the z 2 -type Mn-orbitals in the (ac) plane are shown in the inserts. 
In the 00 phase each Mn-site is statistically occupied by Mn 3+ and Mn 4+ . In the CO-phase the 
fully occupied Mn 3+ -sites are shown with z 2 -type orbitals, while the sites occupied by both Mn 4+ 
and Mn 3+ are indicated by the circles. See the text for the details. 
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FIG. 9: (color online) Bond length as a function of temperature in (Lai_ y Pr y )o.7Cao.3Mn03 for 
the Pr concentration y=1.0, 0.9, 0.85, 0.8, 0.5. The bond lengths for the sample with y=0.7, are 
similar to ones for the y=0.8 and are shown in Fig. 8. 16 and 18 0-samples are shown by closed 
and open symbols respectively. 
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FIG. 10: (color online) (a) The lattice micro-strains as a function of Pr concentration (sp. 
Pnma) at T = 290 K. The anisotropic micro-strain ( JcZ S f / d) [ioo] * s a ^ so shown for T=15 K. 
and 18 0-samples are shown by closed and open symbols respectively. 
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FIG. 11: (color online) Temperature dependences of the integrated intensities of the selected 
magnetic diffraction peaks indicated in the Fig. 4 in the samples with (a) y = 0.5 and (b) y = 0.8. 
The data were collected on heating. 16 O and 18 0-samples are shown by closed and open symbols 
respectively. The lines for the AFM intensities are guides to the eye. For the FM ones the lines 
are fits to the formula described in the text. 
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FIG. 12: (color online) Effective (a) ferro- mp and (b) antiferromagnetic moments tuaf at T=15 K, 
(c) Neel Tn and Curie Tc transition temperatures determined from the ND data as a function of 
y in (Lai_yPr y )o.7Cao.3Mn03. The lines are guides for the eyes. 16 and 18 0-samples are shown 
by closed and open symbols respectively. Neel temperatures for the samples with different oxygen 
mass are the same. 
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FIG. 13: (color online) Effective antiferromagnetic function of the effective ferromagnetic 

moment m 2 F for all the samples with y as a parameter indicated for each point. The top-x and right- 
ly axes show respective phase fractions. The line is a linear fit, y=0.9 and y=1.0 ( 18 0) are excluded 
from the fit. The closed and open symbols represent the 16 and 18 samples, respectively. 
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